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[l]  A  three-dimensional  absolute  velocity  field  on  the 
isopycnal  surface  in  the  northwest  Pacific  subtropical  region 
in  the  summer  is  calculated  from  the  Navy’s  Generalized 
Digital  Environmental  Model  (GDEM)  climatological 
temperature  and  salinity  data  on  a  0.5°  x  0.5°  grid  using 
the  P-vector  method.  GDEM  for  the  area  was  built  from 
historical  (1930-1997)  temperature  and  salinity  profiles.  The 
basic  current  system,  consisting  of  North  Equatorial  Current, 
Kuroshio  and  its  recirculation.  Subtropical  Countercurrent 
(STCC),  is  identified.  STCC  occurs  from  the  surface  to  the 
level  of  CT0  =  25.8  between  18°N  and  23.5°N,  with  a 
maximum  velocity  around  0.1  m/s.  At  ere  =  23.5,  STCC 
‘  originates  east  of  Luzon  Strait  at  around  122.5°E  and 
meanders  eastward.  As  09  increases,  its  origin  shifts  toward 
northeast.  Comparison  between  velocity  and  potential 
vorticity  fields  leads  to  the  conclusion  that  the  STCC  core 
is  above  the  southern  boundary  of  the  subtropical  mode 
water.  INDEX  TERMS:  4576  Oceanography;  Physical:  Western 
boundary  currents;  4572  Oceanography;  Physical;  Upper  ocean 
processes;  4520  Oceanography;  Physical:  Eddies  and  mesoscale 
processes;  4512  Oceanography:  Physical:  Currents;  4532 
Oceanography:  Physical:  General  circulation.  Citation:  Chu, 
P.  C.,  R.  Li,  and  X.  You,  Northwest  Pacific  subtropical  counter- 
current  on  isopycnal  surface  in  summer,  Geophys.  Res.  Lett., 
29(17),  1842,  doi:  10. 1029/2002GL01483 1 ,  2002. 

1.  Introduction 

[2]  One  expects  a  westward  flow  in  the  subtropical 
northwest  Pacific  (Figure  1)  feeding  the  western  boundary 
current  based  on  the  classical  wind-driven  theory.  However, 
early  observations  [e.g.,  Uda,  1955]  show  the  occurrence  of 
eastward  flow  segment  north  of  15°N.  Yoshida  and  Kido- 
koro  [1967a,  1967b]  predicted  the  existence  of  an  eastward 
current  in  the  subtropical  region  especially  20°-25°N  after 
computing  the  Sverdrup  transport  currents  from  the  surface 
wind  stress,  and  named  it  as  the  Subtropical  Countercurrent 
(STCC).  This  current  was  subsequently  confirmed  by  Uda 
and  Hasunuma  [1969]  from  direct  current  meter  observa¬ 
tions  and  geostrophic  calculations,  and  by  White  et  al. 
[1978]  using  historical  expendable  bathythermograph 
(XBT)  data.  Qiu  [1999]  inferred  temporal  and  spatial 
variability  of  STCC  using  the  altimeter  data  from  the  first 
5  1/4  -yr  TOPEX/POSEEDON  mission,  and  found  that  the 
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eddy  kinetic  energy  of  STCC  is  maximum  in  April/May  and 
minimum  in  December/January. 

[3]  Regardless  the  advancement  of  the  knowledge  on 
STCC,  there  is  still  lack  of  quantitative  description  of  STCC 
such  as  its  origin  and  spatial  variability  especially  on  the 
isopycnal  surface.  In  this  study,  the  P-vector  method  on  the 
isopycnal  surfaces  [Chu,  1995,  2000;  Chu  and  Li,  2000]  is 
used  to  calculate  the  absolute  velocity  in  the  northwest 
Pacific  [15°-35°N,  120°-145°E]  from  the  Navy’s  Gener¬ 
alized  Digital  Environmental  Model  (GDEM)  climatologi¬ 
cal  June  temperature  and  salinity  data  set  with  0.5°  x  0.5° 
resolution  [Teague  et  al,  1990].  The  P-vector  method  is  a 
two-step  inverse  method  with  determination  of  the  direction 
of  the  current  (P-vector)  at  first  and  then  the  magnitude 
using  the  thermal  wind  relation.  The  two  necessary  con¬ 
ditions  for  validity  of  inversion  are  easily  cooperated  into 
the  P-vector  method.  The  rest  of  the  paper  is  outlined  as 
follows:  Section  2  depicts  the  Navy’s  GDEM  data  set. 
Section  3  describes  the  computation  in  the  isopycnal  coor¬ 
dinate  system.  Section  4  depicts  the  June  mean  calculated 
absolute  velocity  fields.  Section  5  presents  the  spatial 
variability  of  STCC  in  June.  Section  6  discusses  the  con¬ 
nection  of  STCC  and  the  potential  vorticity  (PV)  field. 
Section  7  presents  the  conclusions. 

2.  The  Navy’s  GDEM 

[4]  Data  for  building  the  current  version  of  GDEM 
climatology  were  obtained  from  the  Navy’s  Master  Oceano¬ 
graphic  Observational  Data  Set  (MOODS),  which  has  more 
than  5  million  profiles  during  1920-1998.  The  basic  design 
concept  of  GDEM  is  the  determination  of  a  set  of  analytical 
curves  that  represent  the  mean  vertical  distributions  of 
temperature  and  salinity  for  grid  cells  (0.5°  x  0.5°)  through 
the  averaging  of  the  coefficients  for  the  curves  found  for 
individual  profiles  [Teague  et  al.,  1990].  Different  families 
of  representative  curves  have  been  chosen  for  shallow,  mid¬ 
depth,  and  deep-depth  ranges,  with  each  chosen  so  that  the 
number  of  parameters  required  to  yield  a  smooth,  mean 
profile  over  the  range  was  minimized.  As  pointed  by  Teague 
et  al.  [1990],  large-scale  oceanographic  features  are  sim¬ 
ilarly  represented  in  both  GDEM  and  the  NOAA  Climato¬ 
logical  Atlas  of  the  World  Ocean  temperature  and  salinity. 

3.  Isopycnal  Coordinate  System 
3.1.  Vertical  Discretization 

[5]  From  the  GDEM  T,  S  data  set,  the  potential  density 
varies  from  22.2  to  27.725.  We  discretize  <t0  with  the 
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Figure  1.  Geography  and  bottom  topography  of  the 
Northwestern  Pacific. 


increment,  Act9  =  0.025  (kg  m~3).  Thus,  we  have  222  de¬ 
layers.  Within  each  layer,  the  density  is  vertically  uniform. 
In  order  to  well  resolve  isopycnal  surfaces,  we  use  cubic 
spline  to  interpolate  the  T,  S  data  into  246  z-levels  with 
three  different  increments:  5  m  from  0-  to  100-m  depth,  10 
m  from  100-  to  1000-m  depths,  20  m  from  1000-  to  2500-m 
depths,  and  50  m  below  2500-m  depths.  Thus,  we  build  up  a 
high-resolution  data  set  [T(zj),  S(zj),  de(z,)]  in  z-coordinate. 
Note  that  the  value  of  27.725  is  the  maximum  value  for  a@, 
computed  from  the  GDEM  T,  S  data  set. 

3.2.  Transformation  of  T,  S  Data  from  z-Coordinates 
to  cre-Coordinate  System 

[6]  The  z-coordinate  potential  density  data  69  (zy)  is 
compared  to  the  discrete  a-values,  a(k).  The  sigma  value 
at  the  bottom  of  the  k-th  isopycnal  layer,  a*(k),  is  used  for 
the  comparison, 

°*(*)  =  5  K*) +  «(*+!)]  0) 

The  depth  for  the  bottom  of  the  a  (k)-layer  is  obtained  by 

4a)  =  Zy,  if  ae(zy)  =  ab(k)  (2) 

and 

4a)=  Zy  +  .  CT/^  7  c  (zy+l  -  Z;)if  ^(zy)  <  cb{k)  <  0e(zy+i) 

at\zj+\ )  —  09{Zj) 

(3) 

where  the  superscript  ‘b’  indicates  the  bottom  of  the  k-th 
isopycnal  layer.  The  thickness  of  the  k-th  isopycnal  layer  is 
obtained  by 

40)  =2*-i  -2*- 

4.  Velocity  Vectors  at  Isopycnal  Surface 

[7]  After  h *(<T)  is  obtained,  we  may  compute  PV  on  the 
isopycnal  surface,  and  then  use  the  P-vector  method  to 


compute  the  absolute  velocity  [Chu  and  Li,  2000;  Chu, 
2000].  It  is  noticed  that  only  the  geostrophic  velocity  is 
obtained  using  the  P-vector  method  [Chu,  1995].  Near 
the  surface  the  geostrophic  current  is  strong  at  least  in 
the  Kuroshio.  Since  the  isopycnal  cr6  =  23.5  outcrops  at 
26°-30°N  (Figure  l)(Figure  2a),  the  isopycnal  surface 
circulations  near  09  =  23.5  might  not  be  consistent  with 
the  direct  measurements  with  current  meters  south  of 
Honsu.  Thus,  the  velocity  at  <j9  =  25.0  is  chosen  to 
represent  the  Kuroshio  Current  and  its  extension  near 
Honsu. 

[8]  Most  important  features  from  the  computation  are  the 
westward  flowing  North  Equatorial  Current  (NEC),  the 
northeast  flowing  Kuroshio  Current  (KC)  and  the  associated 
eddies  and  recirculation  current  from  Kuroshio  (RC),  and  the 
eastward  flowing  STCC.  At  09  =  23.5,  NEC  is  very  weak. 


Figure  2.  Absolute  velocity  vectors  (unit:  cm/s)  in  June  on 
isopycnal  surface  (a)  ct9  =  23.5,  and  (b)a9  =  25.0. 
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Figure  3.  East-west  absolute  velocity  (unit:  cm/s)  in  June 
at  meridional  cross-sections:  (a)  127.5°E,  (b)  129.5°E,  (c) 
133.5°E,  and  (d)  138.5°E.  Here,  the  shaded  part  (negative 
values)  refers  to  the  westward  velocity. 

The  eastward  flowing  STCC  originates  mainly  from  the 
southward  turning  of  RC  (westward  flowing)  at  (127— 
128°E,  23°N).  At  a0  =  25.0  (Figure  2b),  NEC  is  quite  strong. 
Between  128°-134°E  a  weak  current  flowing  eastward 
along  21°N  is  the  north  flank  of  an  anticyclonic  eddy  located 
at  18°-21°N,  128°-134°E  (Figure  2b).  Evident  eastward 
•  flow  (i.e.,  STCC)  occurs  east  of  135°E,  which  is  different 
from  STCC  at  ne  =  23.5.  At  deeper  level  (ct0  =  25.8),  there  is 
no  eastward  current  in  the  subtropics  (figure  not  shown). 

[9]  Two  theories  are  available  to  explain  the  STCC 
formation:  (a)  the  Ekman  convergence  [e.g.,  Roden,  1980] 
induced  by  westerlies  to  the  north  and  trade  winds  to  the 
south,  and  (b)  the  geostrophic  convergence,  which  causes 
the  midgyre  front  and  its  associated  zonal  current  [Cush- 
man-Roisin,  1984;  Kubokawa,  1997],  Since  the  velocity  is 
wholly  determined  by  the  density  structure,  the  identifica¬ 
tion  of  STCC  here  tends  to  confirm  the  validity  of  the 
geostrophic  convergence  theory. 

[10]  Comparison  between  Figures  2a  and  2b  shows  a 
strong  vertical  shear  in  STCC.  The  baroclinic  instability 
associated  with  the  vertical  shear  causes  STCC  to  meander 
between  19°-25°N  with  the  main  axis  at  around  21.5°N, 
and  makes  the  eddy  generation  along  it.  STCC  becomes 
wider  with  a  maximum  width  of  six  degrees  of  latitude 
(18°-24°N)  at  136°-140°E  as  it  flows  eastward  entraining 
water  from  eddies. 

5.  Spatial  Variability  of  STCC 

[11]  The  spatial  variability  of  STCC  is  represented  by 
four  meridional  cross-sections  (127. 5°E,  129. 5°E,  133.5  E, 
138. 5°E)  of  the  east-west  absolute  velocities  on  isopycnal 
surfaces  (Figure  3).  STCC  is  the  eastward  flow  near  the 
Tropics  of  Cancer  (23.5°N),  and  occurs  in  the  upper  layer 
above  o0  =  25.8.  The  core  of  STCC  is  not  located  at  the 
surface.  The  axis  of  STCC  shifts  northward  as  erg  increases, 
and  the  speed  of  STCC  increases  eastward.  For  example, 
the  maximum  velocity  of  STCC  is  around  4  cm  s  1  near 
20°N  between  cj0  —  23.0  and  cr0  =  24.0  at  127.5°E  cross- 
section  (Figure  3a),  greater  than  6  cm  s"1  near  21°N  and 
a0  =  24.2  at  129.5°E  cross-section  (Figure  3b),  greater  than 
9  cm  s_1  near  21.5°N  and  a0  =  24.3  at  133.5°E  cross¬ 


section  (Figure  3c),  and  greater  than  12  cm  s  near  23  N 
and  ct0  =  24.5  at  138.5°E  cross-section  (Figure  3d).  At  the 
ct9  level  where  the  core  of  STCC  located,  the  speed  of 
STCC  is  usually  higher  than  the  speed  of  NEC.  This  is 
because  the  core  of  NEC  is  located  at  deeper  layers.  The 
volume  transport  of  STCC  at  138. 5°E  is  estimated  as  8  Sv 
(1  Sv  =  106  m3  s_1). 

6.  Potential  Vorticity 

[12]  STCC  mainly  comes  from  southward  turning  of  the 
recirculation  current  from  Kuroshio  (i.e.  RC),  which  orig¬ 
inates  at  southern  latitudes  (~15°N)  where  PV  is  relatively 
low.  Scale  analysis  shows  that  the  dissipative  force  alone  is 
not  sufficient  to  remove  the  low  PV  anomalies  [Pedlosky, 
1987]  for  a  narrow  boundary  current  such  as  KC  and  RC. 
Talley  [1988]  pointed  out  the  existence  of  a  minimum  PV 
area  corresponding  to  the  subtropical  mode  water  located  at 
(je  =  25.4  level  in  the  western  Pacific.  The  minimum  PV 
enclosed  by  25  x  10_1°  cm_1s_1  at  the  cross-section  of 
138.5°E  (Figure  4b)  indicates  the  location  of  the  subtrop¬ 
ical  mode  water.  The  STCC  core  is  located  above  the 
southern  boundary  of  the  subtropical  mode  water  (i.e., 
minimum  PV)  from  the  comparison  between  Figures  3d 
and  4.  This  result  consists  with  the  recent  modeling  studies 
[e.g.,  Kubokawa,  1997]. 

7.  Conclusions 

[13]  The  goal  of  this  study  is  to  investigate  the  spatial 
variability  of  summer  Subtropical  Countercurrent  on  iso¬ 
pycnal  surfaces  in  the  northwestern  Pacific  using  die  Navy’s 
unclassified  GDEM  data  on  0.5°  x  0.5°  resolution.  We  use 
the  P-vector  inverse  method  for  determining  the  absolute 
geostrophic  velocity  on  the  isopycnal  surface.  The  results 
are  summarized  as  follows. 

1.  The  eastward  currents  at  upper  levels  (above  cr0  = 
24.5)  reveal  multi-core  structure.  STCC  is  located  at  19°- 


Figure  4.  Meridional  cross-section  (138. 5°E)  of  June 
potential  vorticity  (in  10“ 10  erne's-1). 
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23  °N.  Eastward  flow  segments  occurring  south  of  17°N  are 
usually  the  southern  flanks  of  local  cyclonic  eddies  rather 
than  a  part  of  a  continuous  current.  Below  the  level  of  ere  = 
24.5,  STCC  has  a  single  core  located  19°-23°N  west  of 
134°E  and  22°-26°N  east  of  134°E. 

2.  STCC  originates  at  122.5°E  at  ere  =  23.5.  As  ag 
increases,  the  origin  of  STCC  shifts  eastward.  The  major 
source  of  STCC  is  the  southeastward  turning  of  the 
recirculation  current  from  Kuroshio.  STCC  strengthens 
and  broadens  as  it  flows  eastward.  The  volume  transport  of 
STCC  at  138.5°E  is  estimated  as  8  Sv  (1  Sv  =  106  m3  s-1). 

3.  The  STCC  core  is  not  located  at  the  surface,  but  the 
layer  between  eg  —  23.0  and  erg  =  24.7,  where  STCC  is 
stronger  than  NEC.  Moreover,  STCC  core  is  above  the 
southern  boundary  of  the  subtropical  mode  water. 

4.  STCC  disappears  below  the  level  of  ag  =  25.8,  that  is 
to  say,  STCC  is  an  upper  ocean  current.  Its  maximum  depth 
is  around  300  m. 
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